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Abstract 

The spectral properties of the Probenius-Perron operator of one-dimensional 
maps are studied when approaching a weakly intermittent situation. Nu- 
merical investigation of a particular family of maps shows that the spectrum 
becomes extremely dense and the eigenfunctions become concentrated in the 
vicinity of the intermittent fixed point. Analytical considerations generalize 
the results to a broader class of maps near and at weak intermittency and 
show that one branch of the map is dominant in determination of the spec- 
trum. Explicit approximate expressions are derived for both the eigenvalues 
and the eigenfunctions and are compared with the numerical results. 
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I. INTRODUCTION 



Correlation functions play an important role in the characterization of chaotic systems. 
In certain classes of systems they decay exponentially. Cases when it is proven are Axiom A 
systems [^]-|^, mixing maps ^ and maps which are analytic on a set of rectangles 0. It is 
generally believed that hyperbolic systems usually exhibit an exponential correlation decay, 
although the rigorous proof of this can be quite a difficult task in particular cases . On the 
other hand there are similarly important classes of systems where the decay of correlations 
is slower than exponential. Such a behavior has been found in intermittent one-dimensional 
maps and in Hamiltonian systems with mixed phase space f^. Intermittency, the alter- 
nation of almost regular and chaotic motion, can precede birth of a pair of stable-unstable 
fixed point or periodic orbit (tangent bifurcation) P,p^, can be caused by a marginally 
unstable fixed point |1TT|,|7|JT^-[14| or the marginal behavior at the boundary of an island in 
general Hamiltonian systems [^. It can also arise when a chaotic attractor becomes unsta- 



ble due to crises [jT5[, effect of noise or driving of another attractor |[16|JT7[] . Intermittency is 
known to possess several special properties. The spectrum of Renyi-entropies shows a phase- 
transition-like behavior |[T3| jT^ -pT| , the decay properties of truncated entropies, complexity 
and the dynamical fiuctuations also differ from those of expanding maps p!^ , |20|p^ -p5[ . The 
kind of intermittency we study is due to a marginally unstable fixed point |TT|,|^JT^-|T3|. The 
vicinity of the fixed point is visited by typical trajectories and the time needed to leave this 
vicinity scales with a power of the distance from the fixed point instead of the logarithm as 
for an unstable fixed point. Depending on the actual map there can even exist a smooth 
invariant density to which initial densities converge, as in case of our example. This case is 
called weak intermittency, otherwise we talk about strong intermittency 0. 

The decay of correlations may be related to the spectral properties of the Frobenius- 
Perron operator describing the time evolution of the probability density given in the phase 
space of the system in question. The rough picture is that the largest eigenvalue is unity 
(except for transient chaos |^ , not to be treated here) and the corresponding eigenfunction 
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is related to the natural measure, while the leading asymptotics of the correlation decay is 
related to a spectral gap between the largest eigenvalue and the rest of the spectrum. Note 
that in two or higher dimensional systems even in hyperbolic cases the direct study of the 
spectrum of the Frobenius-Perron operator involves a series of difficulties, related to both 
the suitable definition of the Frobenius-Perron operator and the relevant function space in 
which the eigenvalue problem is to be solved. 

In case of one-dimensional maps the Frobenius-Perron operator takes the form Lip{x) = 
Y,z:f{z)=x ■ Here the sum goes over the preimages z of x. Existence of a unique invariant 
density (belonging to eigenvalue Aq = 1) is proven for expanding piecewise maps |^ and 
fully developed maps with negative Schwarzian derivative |2^. If a unique ergodic measure 
exists the correlation function can be expressed in terms of the Frobenius-Perron operator. 
If the spectrum is discrete eigenvalues A„ and eigenf unctions ipn may be used to explore the 
time evolution of an initial probability density (p{0,x). Then decay of correlations and the 
convergence of (p(t,x) towards the invariant density is governed by Ai (where Aq > |Ai| > 
IA2I > ... > |A„| >...). 

The relation between the spectral properties and the correlation decay becomes much 
more complex if the latter one is slower than exponential. In such a case, arising typically in 
intermittent systems, the largest eigenvalue must be an accumulation point of the spectrum. 
The standard methods developed for determining the eigenvalues of operators 



similar to the Frobenius-Perron operator can not be used to determine the non-leading 
eigenvalues when the system is close to an intermittent situation. A direct study of the 
spectrum of generalized Frobenius-Perron operators in case of intermittency can be found 
in Ref. |Q. Working in the space of functions with bounded variation a nonzero essential 



spectral radius ||3^ is obtained, which reaches Aq for the case of L . Study of the spectrum 



for intermittency that precedes tangent bifurcation is presented in |^ . 

Our motivation has been to calculate the spectral properties near and at intermittency 
and to establish their link to the correlation decay. In the present paper we concentrate 
on the first task, while the derivation of the correlation functions will be published else- 



where |3^. In order to avoid the difficulties arising in higher dimensional systems we con- 
sider piecewise analytic one-dimensional maps and choose the function space of real analytic 
functions. Note that other methods and other function space has been used than in Ref. |p3|| , 



which allows us a more detailed description. The paper is organized as follows. In Section 
2 we study the spectral properties of the Frobenius- Perron operator numerically in a family 
of fully developed chaotic maps |37|j38|| , where the parameter is tuned up to weak intermit- 



tency. Some aspects of the numerical work are discussed in Appendix A. Section 3 contains 
the analytical results, which are valid for a wider class of maps. Details of the calculation 
are given in Appendix B. 

II. A MODEL AND NUMERICAL CALCULATION 



Our model is the piecewise parabolic map [38|. This is probably the simplest dynamical 



system which can show intermittency. The general form of the map is 



1 + ^ - \/ 1 - ^ + 4r 1 - 2x\ 

fr{x) = ^ . (1) 

It maps the interval [0, 1] twice onto itself for every value of the parameter r in the interval 
[—1, 1] . Its name comes from the fact that its inverse has parabolic branches. The point 
X = is a fixed point which becomes marginal in the intermittent case r = 1 p8i|39| ,[7|. 



For the other special case r = we obtain the uniformly expanding tent map. A special 
property of the piecewise parabolic map is that the stationary probability density has a 



simple form ||3^, Pr{x) = 1 + r(l — 2x). This is still normalizable for r = 1, which means 
that the map possesses weak intermittency. Note that the average Ljapunov exponent is 
positive. 

In order to study the eigenvalue spectrum of the Frobenius-Perron operator we derive 
a matrix representation on the basis of the Legendre polynomials P^. After normalization 



Qkix) = yk+l/2Pk{x) these polynomials form a complete ortho normal system on the 
interval [— 1, 1] . In order to use them we have to transform the map and this basis to the 
same interval. 



In the case of our model the action of the Probenius- Perron operator on an arbitrary 
function (p{x) can be given by 

L<p{x) = [^{F,{x)) + ^{F,{x))]\F^{x)\ , (2) 

where -^1,2(2^) — | T (| — + |a;^^ stands for the inverse branches of the map. 

Expanding the function (p{x) in terms of Qk{x) as (p{x) ~ J2kLoi>kQk{x) we obtain the 
(infinite) matrix representation of the Frobenius-Perron operator, 

CO 

Lip{x) = J2 LikbkQi{x) , (3) 

k,l=0 

where the matrix elements L;^ are of the form Li^ — J^i Qi{x)LQk{x) dx . In the numerical 
calculation we have truncated this matrix confining the indices between and N . The 
integrals have been evaluated by numerical integration. Exploiting the symmetry of the 
map the matrix size has been effectively reduced by a factor of 2. Then the eigenvalue 
spectrum A^^\ n = 0, 1, 2, ... of the matrix has been determined. It was quite remarkable 
that the eigenvalues turned out to be always positive real numbers. This calculation was 
performed for different values of the parameters r and N . Fig. 1 shows the spectrum as a 
function of the parameter r at two different values of N. The solid parts of the curves in 
Fig. 1 correspond to cases where, at the given value of iV, the precision of the eigenvalues 
have reached 0.003 (which is approximately the resolution of the figure). The precision was 
determined by comparing results for different values of N . 

It is seen that all the eigenvalues A„ with fixed n tend to value 1 for r — >■ 1. On the other 
hand we find that near any fixed A value the density of the eigenvalues increases suggesting 
a continuous spectrum at intermittency. The eigenfunctions of the operator also show a 
remarkable behavior when approaching the intermittent case r = 1. Apart from the first 
one, they become more and more concentrated in the vicinity of the fixed point (see Fig. 2). 
That means, the limiting eigenfunctions are singular. One might expect the strange case 
that the eigenvalue Ao = 1 becomes infinitely degenerate in the space of analytic functions. 
This is avoided since the eigenfunctions become singular in the limit. At the same time the 
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tails of the eigenfunctions outside the vicinity of the fixed point are smooth functions and 
become more and more similar to the invariant density p{x) (see inset in Fig. 2). In fact, 
with a suitable normalization they converge to p{x). 

It has been observed that for r < 0.9 the eigenvalues converge quite fast when increasing 
the value of N. Going closer to the case r — 1 the convergence becomes slower and the 
limiting eigenvalues get closer to the largest one (Aq = !)• In the case r = 1 the eigenvalues 
converge to unity, but slower than exponentially. At the same time the eigenvalues with 
higher index converge slower, consequently the closer r is to unity the less eigenvalues are 
precise. Its reason is twofold. On the one hand, the spectrum becomes denser and denser 
when approaching the intermittent case, which certainly leads to numerical problems. On 
the other hand, the eigenfunctions (pni^) close to the intermittent situation after suitable 
conjugation can be well approximated by powers x"'~^. As it is explained in Appendix A 
this leads to the problem, that only a few eigenvalues and eigenfunctions can be obtained 
numerically. 

III. ANALYTICAL RESULTS 

We can understand the behavior of the spectrum near the intermittent case by keeping 
only the term of the Probenius- Perron operator L(p{x) = J2z:f{z)=x ^^^^ corresponds 
to the lower inverse branch. This branch contains the intermittent fixed point. (We shall 
denote the corresponding operator by Li, while the operator containing the upper inverse 
branch of the map will be denoted by L2.) The first explanation is the following. The 
numerical results show that the eigenfunctions, except for the invariant density p{x), are 
concentrated in the vicinity of the fixed point. Therefore the most important region is this 
neighbourhood. The contribution of L2 comes from the rightmost piece of the interval [0,1], 
while this piece is reached from the neighbourhood of the maximum in one iteration. At the 
same time the eigenfunctions are smooth and (except for p{x)) have quite low value there 
compared to the vicinity of the fixed point. Therefore the contribution of the upper branch 
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is small and is a slowly varying function. So we first study Li and later we return to the full 
Frobenius-Perron operator to estimate the effect of the upper branch. The considerations are 
quite general, we only assume that the map is single-humped, we assume special properties 
of the map at the endpoints and the maximum, and assume that the map is smooth and 
expanding everywhere else. 

Without restricting generality we let the map act on the interval [0, 1] so that the fixed 
point is at X = 0. We assume that the map satisfies the following formulas such as model 
(0) does: 

fix) = (1 + e)x + ax^ + 0{x^) if x < 1 , (4) 
f\x) = d{l -x) + - xf) if 1 - X < 1 , (5) 
Fi,2(x) =XT (^i,2£(l -x) + /ii,2(l - xf) 

+ 0{{l-xf) if 1 - X < 1 , (6) 

where x is the location of the maximum, e,a,d,gi^2,hi^2 are suitable positive constants 
and Fi^2 are the inverse branches of /(x). For e —>■ keeping a,d,gi^2,hi,2 fixed we are 
approaching the intermittent situation. The eigenvalue equation for Li takes the form 

A^(/(^)) = ^, 0<x<x. (7) 

We can solve it in the neighbourhood of the origin by introducing uj{x) = \n(f{x) . Since 
/(x) — X is small when we are close to the intermittent case uj{f{x)) — u!{x) can be approx- 
imated by the linear term of Taylor expansion. Then one obtains 

Substitution of @ and integration yields 

where i/(A) = —j^ — 1 + |. Let us restrict ourselves to the case of analytic eigenfunctions, 
that means i^(A) = n — 1, where n is a positive integer. Then 



The corresponding eigenvalue is obtained from the expression below (H) as A„ ~ e""*^. 

It is interesting to see that even the exact eigenvalues of Li (see(^) for a not necessarily 
small e can be obtained from the stability of the fixed point. The reason is the following. If 
we know the eigenfunction v?(x) in any small neighbourhood Vq of the fixed point, then by 
we can determine it in the image Vi of this interval. Repeating this procedure, say k times, 
the part of the interval Vk lying outside the previous one Vk-i will be mapped outside . 
Hence we can find the eigenfunction satisfying (|^) in [0, x]. For small enough neighbourhood 
Vo the function f{x) can be approximated linearly, which yields eigenfunctions x""^ and 
eigenvalues 

An = (1 + 5)-", n = l,2,... . (11) 

Moreover, these are the exact eigenvalues of Li, since the relative precision of the linear 
approximation can be arbitrarily increased by shrinking the interval Vq. However, for the 
determination of the eigenfunctions the quadratic term in (^ is important. 

After these considerations we can estimate the effect of the upper inverse branch of the 
map. The eigenvalue equation may be written as 

\f[{xMf{x)) = ip{x)+x{x), (12) 
x{x) = ^{F,U-{x)))f[ix)F^{f{x)). (13) 

It is seen from Eq. ([To|) that very close to the fixed point the eigenfunction starts as 
fn{x) ~ f3x^~^ (we do not fix the normalization here, f3 may contain a power of e) and 
it is concentrated in a region of width of order e. We assume that the eigenfunction gets 
small corrections compared to Eq. (|10D- After estimating the value of ipn and writing down 
power-series forms for it in the neighbourhood of a; = and x we obtain results for A„ and 
coefficients of the expansion V9„(x) = I]fclo A^'^ • 

(3^ = C(/3e"-'") , if m < n - 1 , (14) 

A„ = (l+£)-'^ + 0(£2). (15) 
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The actual value of the correction depends on the shape of the map but its order of magnitude 
does not. The order of the correction in (^) was also shown by the numerical results for 
the map (|1]) for the n = 1, 2, 3 eigenvalues. For the perturbation of the eigenfunction (see 
(p^)) the case m = shows that the eigenfunction is shifted, and the ratio of this shift (3o 
and the maximum is of order e. This was numerically verified for n = 2, 3 (see the nonzero 
starting value in Fig. 3). The value (3i gives the initial slope of the eigenfunction at the 
fixed point. j3i can be compared to the slope C(/3e"~^) of the straight line connecting the 
origin and the maximum of Their ratio is also proportional to e, what was numerically 
checked for n = ?>. 

It is important to mention, that not only values of A„,/3o and [3i show good agreement 
between the analytic and numeric results. The shape of the eigenfunction itself agrees very 
well, too. The agreement is surprisingly good for r = 0.999 . Fig. 3 also supports the above 
estimations, namely, that the deviation, which is due to the effect of the upper branch, is 
proportional to e. One should notice, that we did not get the eigenvalue Aq = 1 this way. 
The reason is, that for n = the eigenfunction is not concentrated in the vicinity of the 
fixed point, as indicated in the beginning of this section. This means, the contribution x{^) 
of the upper inverse branch of the map is not negligible in this case. 

In the above estimation of the effect of the second branch n was assumed to be constant. 
However, the effect of L2 remains also small in the limit r ^ 1 when n tends to infinity in 
such a way that A„ is kept constant. This can be realized using a series of r„ values which 
satisfy A„(r„) = A. The effect of L2, i.e. the difference 5\n of analytic eigenvalues ( pT]) for 
Li and numeric eigenvalues for L was determined for several values of n and r. Since the 
spacing AA„ = A„ — A„+i of the eigenvalues decreases in the limit the relative precision 
5A„/AA„ is important. Fig. 4 shows that this ratio is roughly constant in function of n 
(depends on A only) when it is considered as a function of n and A, and it is small for A ~ 1. 
According to these numerical findings one may assume that also the expression (p^ ) for the 
eigenf unctions is valid with a fixed precision in the limit e — » for fixed A. Therefore an 
approximating eigenfunction can be obtained for the intermittent case r = 1 as the limit 
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of ([lOD, which can also be calculated using (|]), the expression of z/(A) and keeping A fixed. 
This way one obtains 

^x) ^ A,-^e-'-^^'^^ , (16) 
while the spectrum is found to be continuous. 
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APPENDIX A: REPRESENTING THE EIGENFUNCTIONS ON AN 

ORTHONORMAL BASIS 

In this appendix we explain why an orthonormal basis is not well suited to get many 
eigenvalues correctly near an intermittent state. This is easier to see after applying the 
following conjugation to the piecewise parabolic map: 

y = h{x) = ^ (i + -) ■ (Ai) 

X + ^ \ a) 

Then the eigenfunctions of the map can be well approximated by powers near intermit- 
tency. Therefore (<l>m|-^|y") is a good matrix representation of L, where {$n} is the adjoint 
basis to {?/"}. 
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Using an orthonormal basis {qn{y)} the matrix elements Lnm = {(ln\L\(lm) can be ob- 
tained by the following transformation: 

{^,nW)=T.W-lLklWir^, (A2) 
l,k 

where Wnm = Jo Qniy)y"^ dy . How many eigenvalues we can get correctly depends on how 
weakly conditioned the matrix Wnm is- 

In our case {qn{x) is the normalized Legendre polynomial transformed to [0, 1]) the 
matrix elements can be calculated analytically using Rodrigues formula and applying partial 
integration: 



y/2n + l(m! 



|2 

for n < m 



= <! ("^~"')'("^ + ^ + (A3) 
for n > m 

The maximum number Nmax of eigenvalues and eigenvectors that can be reliably computed 

can be approximately read off from the singular value decomposition representing a matrix 

W as Wij = UiiDiiV^y Here U and V are orthogonal n x n matrices and D is diagonal one 

with the property Da > Djj > for i < j. We have determined Nmax from the condition 

Dm at > 6, Dm 117V +1 < 5 (A4) 

with two different choices for 6 (10~^ and 10~^°). It is seen in Table 1, that Nmax increases 
slowly with increasing n. Therefore we conclude that close to intermittency the shifted 
Legendre polynomials represent a good basis to determine the first few eigenvalues and 
eigenfunctions. Beyond that they are expected to produce misleading results. We surmise 
that these conclusions are intrinsic for any orthogonal basis since the matrix W is nearly 
degenerate. This corresponds to the fact that the basis functions are nearly parallel for 
large n when angles are defined through the usual C2 norm. 

APPENDIX B: EFFECT OF THE UPPER INVERSE BRANCH 

To obtain estimations ([T^), (|TB|) we start from Equations (|T2|)-(^). It follows, that 
the order of the maximum value of '^{x) is /^e""^. The numerical calculation not only 
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supports these estimates but shows that the function outside that region is smooth and 
slowly varying. The integral of the eigenfunction should be zero since it should be orthogonal 
to the first left eigenfunction of L. At the same time its integral in the above mentioned 
region near is of order /Je", hence it should be of the same order around the position x 
where / has its maximum. Thus the Taylor expansion ip{x) = J2T=o^k{x — x)*^ possesses 
coefficients of order /^e" at most, Ck = Going to the iterate of x we apply equation 

(|T^) for X ~ f{x) = 1. Using Eq. we obtain estimation for the Taylor series around 
the point x = 1. The second term of (|l^) connects it to (f{x) around x = 0. Since 
the expressions F2{f{x)) and /(x) in (|13]) start linearly for x ~ the coefficients in the 
expansions xi^) = J2kLo CkX^ and ^p{x) = J2kLo Ck{^—x)^ are of the same order of magnitude: 
Co, do = Cfc, 4 = Oi/3e^) , > 0. 

Returning to the neighbourhood of the intermittent fixed point we apply (|12D for x ~ 0. 
We describe the corrected eigenfunction by the expansion ip{x) = Y^kLoPkX^, where Pk, k < 
n — 1 are small and ~ /3, which ensures that the eigenfunction gets small correction. 
Substituting the expansions for Lp{x)^x{^) ^"^^ ® i^ito expanding the powers and 

collecting terms proportional to x™ we obtain 

fc=/nt(f ) 

m-2 

+ [A(1 + - l] /5„ = Cm+Y. ^(/^fc) ' (Bl) 

A:=0 

where hm,k = {{Jlk) + '^{m-k-i})^ (g) = if g > p or g < . Using this equation and 
the estimations for Ck it can be shown recursively, that (3^, = 0{(3e'^~'^) , if m < n — 1 . 
For m = only the term containing (3^ remains on the left hand side. Since the eigenvalue 
(p!T|) is expected to get small correction the coefficient of (3^ in ( piP is of order e . This 
yields /3o = 0{(3e'^). For m > the sum in (pT|) is not suppressed by Cm, so we get 
/3m = Oi^e"^ Pm-i), which proves our statement for (3m- However, for m = n — 1 the 
eigenvalue (|1T]) would give zero for the coefficient of /3m, therefore this coefficient depends 
strongly on the perturbation of A . At the same time in the case m = n — 1 the order of 
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Pm — = /3 is known, hence this case determines the corrected value of A with the help 
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FIGURES 

FIG. 1. Eigenvalues of the Probenius- Perron operator L for the piecewise parabohc map with 

maximal matrix index (a) N = 40, (b) = 80. Solid lines: precision higher than 0.003; dashed 
lines: precision lower than 0.003. 

FIG. 2. The eigenfunction y?2 for r = 0.9 (solid line), r = 0.97 (dashed line) and r = 0.99 
(dotted line). Inset: The eigenfunctions (pQ (solid line), (fi (long dashes), (p2 (short dashes), (p^ 
(dotted line) and (^4 (dashed-dotted line), for r = 0.999 with different normalization. 

FIG. 3. Comparison of eigenfunctions (pi,Lp2,(pz and (/34 of Li obtained analytically (solid lines) 
with corresponding eigenfunctions of L calculated numerically for r = 0.99 (dotted lines) and for 
r = 0.999 (dashed lines) in the vicinity of the fixed point. For identification of the curves we note 
that ifi < < 'Ps < ¥'4 at the right edge of the figure. 

FIG. 4. Effect of the second branch on the eigenvalues compared to the eigenvalue spacing in 
function of A and n. 0:n = l, + : n = 2, □:n = 3, x : n = 4, A:n = 5, *:n = 6. 
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TABLES 



n 


() = !(]-■' 


d = 10-1" 


5 


5 


5 


10 


8 


10 


20 


10 


16 


50 


12 


21 


100 


14 


24 


500 


18 


32 


1000 


20 


35 



TABLE I. Values of Nmax in function of n for different S. 
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0.4 - 



0.3 - 



0.2 - 



0.1 - 




0.75 




0.8 



0.85 



0.9 



0.95 



